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ABSTRACT 

In a planetary system with two or more well-spaced, eccentric, inclined planets, secular interactions 
among these planets may lead to chaos. The innermost planet may gradually become increasingly 
eccentric and/or inclined, as a result of the secular degrees of freedom drifting towards equipartition of 
AMD (angular momentum deficit). This "secular chaos" is known to be responsible for the eventual 
destabilization of Mercury in our own Solar System. Here we focus on systems with three giant 
planets. We characterize the secular chaos and demonstrate th e criterion for it to oc cur, but leave a 
detailed understanding of secular chaos to a companion paper ( |Lithwick fc mil[20T0| . 



After an extended period of eccentricity diffusion, the inner planet 's pericentre can approach the star 
to within a few stellar radii. Strong tidal interactions and ensuing tidal dissipation extracts orbital 
energy from the planet and pulls it inward, creating a hot Jupiter. In contrast to other proposed 
channels for the production of hot Jupiters, such a scenario (which we term "secular migration") 
provides an explanation for a range of observations: the pile-up of hot Jupiters at 3-day orbital 
periods, the fact that hot Jupiters are in general less massive than other RV planets, that they may 
have misaligned inclinations with respect to stellar spin, and that they have few easily detectable 
companions (but may have giant companions in distant orbits). Secular migration can also explain 
close-in planets as low in mass as Neptune; and an aborted secular migration can explain the "warm 
Jupiters" at intermediate distances. In addition, the frequency of hot Jupiters formed via secular 
migration increases with stellar age. We further suggest that secular chaos may be responsible for 
the observed eccentricities of giant planets at larger distances, and that these planets could exhibit 
significant spin-orbit misalignment. 



1. INTRODUCTION 

While around 10% of sun-like stars surveyed harbor 
Jovian-mass planets, only ~ 1% are orbited by so-called 
hot-Jupiter s with periods short - ward of ^ 10 da ys (see 



reviews by Marcy et al. 



Udry fc Santos 2007) 



|2005| 

There appears to be a pile-up of hot Jupiters around 
3 day orbital periods. This excess is genuine and has 
been confirmed by both r adial velocity and transit sur 



veys (Gaudi et al. 2005| 



i05| IButler et^, 
et al.|^n05t|l''ressin et"anp007p . Outward of hot Jupiters 
there appears a deficit of gas g iants with p eriods o f 10 to 



2006 Gumming 



100 days (the "p eriod valley;" |Udry et al . 2003; W itten- 
myer et al.||2010D . 

According to current theories of planet formation, hot 
Jupiters could not have formed in situ, given the large 
stellar tidal field, high gas temperature, and low disk 
mass to be found so close to the star. Instead, hot 
Jupiters most likely formed beyond a few AU and then 
migrated inward. Candidate migration scenarios that 
have been proposed include protoplanetary disks, Kozai 
migration by binary or planetary companions, and scat- 
tering with other planets in the system. While each 
of these mechanisms may have contributed to the hot 
Jupiter population to some degree, the question remains 
as to which is the dominant one. The dominant mech- 
anism has to explain a variety of observed correlations. 
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In Sj2j we review some of these correlations and provide 
a critical assessment of the above three mechanisms. 

In this work, we propose a fourth channel for producing 
hot Jupiters, namely planet migration by secular chaos. 
Secular chaos may arise in planetary systems that are 
well-spaced and are dominated by long-range secular in- 
teractions. A system of two non-coplanar planets can be 
chaotic, but only i f their ecc entricities and inclinations 
are of order unity ( L ibert fc: Hcnrard 2005: Migaszewski 
& Gozdziewski'2009T^ aoz ct al. 2010). So in this contri- 



bution we focus on systems with three planets. The crite- 
rion for secular chaos is less stringent, and the character 
of secular chaos is more diffusive, differing from that of 
the two planet case. This diffusive type of secular chaos 
promotes energy equipartition between different secular 
degrees of freedom. The physics behind s ecular chaos 
is analyze d in detail in a companion paper (Lithwick &: 



Wu 20101, where we show that Mercury, the mnermost 
planet m our Solar system, experiences a similar type of 
secular chaos. Mercury may consequent ly be removed 
from the Solar system (| Laskar 2008; Batygin fc Laughlin 
'2008! |Laskar fc Gast incau 2009). 

Secular chaos tends to removes angular momentum in 
the inner most planet gradually, causing its pericenter 
to approach the star. Tidal dissipation may then re- 
move orbital energy from this planet, turning it into a 
hot Jupiter. Hot Saturns or hot Neptunes may also be 
produced similarly. Such a migration mechanism, which 
we term "secular migration," can reproduce a range of 
observations. It also predicts that in systems with hot 
planets, there are other giant planets roaming at larger 
distances. 
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2. HOT JUPITERS: OBSERVATIONS AND THEORIES 
2.1. Observations 

There is a sharp inner cut-off to the 3-day pile-up of 
hot Jupiters. They appe ar to avoid the regio n inward of 
twice the Roche radius (|Ford &: Rasio||2006p, where the 
Roche radius is the distance within which a planet would 
be tidally shredded. New data spanning two orders of 
magnitude in planetary masses (and including planet ra- 
dius measurements) have strengthened this claim. There 
are only 5 known exceptions lying inward of twice the 
Roche radius, and the rest mostly lie between twice and 
four times the Roche radius. 

Hot Jupiter s appear to be less mass i ve than more dis 



19861. It asserts that the viscous protoplanetary disk 



tant planets (Patzold & Rauer 2002 Zucker & Mazeh 
20021. For planets discovered with the radial velocity 
method, close-in planets have projected masses (Msini) 
less than twice Jupiter's mass, excluding planets in mul- 
tiple star systems. But numerous f urther out planets 
have M sini > 2Mj ( |Udry fc San"tosl|2007, Fig. 5). 

Many hot Jupiters have orbits tnat are misaligned 
with the spin of their host stars. The angle between 
the orbit normal of a transiting planet and the spin axis 
of its host star (the stellar obliquity) can be probed with 



the R ossiter-McLaughlin (R-M) effect (e.g. Winn et al 



2005 1 . Planets were presumably born in a disk aligned 



with the stellar spin. Therefore measurements of the 
present stellar obliquity provide a stringent constraint 
on the migration scenario. Analysis of the first 11 sys- 
tems with R-M measurements found that the majority 
were consistent with perfect ali gnment, while a small mi - 
nority were highly misaligned ( jFabrycky fc Winn||2009 ) . 
However later analysis that included more systems (26 
in total) found that most are misaligne d, and many o f 
these are even in retrograde orbits (Triaud et al. 2010). 
The reason for this discrepancy is currently unclear. 

Hot Jupiters also tend to be alone, at least out to a 
few AU. From radial velocity surveys, ^ 30% of planets 
are in multiple pla net systems (includ ing ones with ra- 
dial velocity trends, Butler et al.||2006 ), while only 5 hot 
Jupiters are (HD 18712 3b, HlPTteTUb , ups A ndb, HAT- 
P-13b and HD 217107b; |Wright et al .l2009; He brard et al. 
|2010a|); i.e. fewer than 10% of hot Jupiters are known 



to have companions within a couple AU. This relative 
deficit also shows up in the transit sample, where most 
attempts at detecti ng transit timing v ariations ca used by 



close companions (Holman & Murr^ 2005; Agol et al 
2005) have been unsuccessful 



El 



Eabus et aE'^DDS 



Csizmadia et al.|2010[ Hru dkova et al.|2010 |, excepFtor , 



perhaps, IVlaciejewski et al. '( 2010b|a ); l-'ukui et al. (2010). 



2.2. Migration Theories 

The successful migration scenario has to explain these 
and other observed correlations. There are two categories 
of migration scenarios. One is to migrate the planet 
within a gaseous disk (disk migration). The other is to 
generate a high eccentricity in the planet, bringing it suf- 
ficiently close to the star that tidal dissipation circular- 
izes and shrinks its orbit into that of a hot Jupiter. The 
latter category includes Kozai migration, planet-planet 
scattering, as well as the secular chaos that we propose 
in this work. 

We first examine disk migration, a theory that pre- 
dated the discovery of hot Jupiters (Lin & Papaloizou 



carrie s the planet inward (Chambers 2009 Rice et al 



20081. The presence of mean-motion resonance pairs 



among observed planets seems to support this scenario. 
However, to produce the observed pile-up of hot Jupiters 
at around 3 day orbital periods, the inward migration 
has to be halted. 



Lin 



As discussed by 
this could be achieved if the disks are truncated 
stellar magnetosphere at a radius that corresponds to 
twice the planet orbital period. Disks are likely trun- 
cated at the co-rotation radius, where the orbital period 
of the disk material equals the spin period of the star. 
The observed distribution of rotation p eriods of pre-main 



seque nce stars appear to be bimodal (Herbst & Mundt 



2005 1 with location of the long period peaks varying trom 
4 to 8 days for different clusters. This could be con- 
sistent with the period distribution of hot Jupiters. It 
is also plausible that hot Jupiters thus migrated would 
tend to be lower in mass, although that has yet to be 
shown. A difficulty with this scenario is that it should 
produce planets whose orbit normals are aligned with the 



stellar spin axis (for an opposite view, however, see Lai 
et al.||2010). In addition, it remains to be argued why 



hot Jupiters rarely have close companions (within a few 
AU) if they are migrated inward by powerful disks. And 
there is no natural explanation in this scenario for the 
avoidance of twice the Roche radius. 

We now turn our attention to planet-p lanet scattering. 
First proposed by Rasio & Ford| (19961, it asserts that 
close encounters between planets can induce extreme ec 



centricity in one of them jF ord et al. 2001 Papaloizou 
' Wask 



k. Terquem 2001 Ford & RasiO||2008j), wHich may then 



. . . iLf 

be tidally circularized to form a hot Jupiter. Such a the 
ory reproduces the observed eccentricity distribut ion of 
(non-hot) e_xtra-solar planets tha t have e > 0.2 ( [Chat 



terjee et al.|2008 Juric & Tremainc 2008 ) , as well as the 



obser ved close-packed nature of pairs (Raymond et al 



20091. It may also account for the high inclinations of 
J 



hot Jupiters. However, it is unclear if the initial condi- 
tion of a compact and highly unstable planetary system 
as required by this theory is applicable to planets emerg- 
ing from a gaseous disk (|Matsumura et al.,,2010) ). Also, 



Chatterjee et al. (2008) find that the inner planets tend 
to be the most massive ones, contrary to the observed 
correlations. Furthermore, since scatterings are sudden, 
it is diffi cult for tides, a slow p rocess, to halt scattered 
planets (Nagasawa et al. 2008). This theory also pre- 
dicts readily detectable outer planets that are responsible 
for scattering and producing the observed hot Jupiters. 
They are, however, not observed. 

Lastly, we comment on Koz ai migration. First pro- 
posed by|Wu & Murray (2003), it asserts that a highly 
inclined co mpanio n star can induce Kozai osci llations 
( |Kozai|[T962, Eggleton fc Kiseleva-Eggleton||200l[ ) in the 
planet, gradually exciting the planet to a high enough ec- 
centricity that it approaches the central star, whereupon 
tidal dissipation circularizes it into a hot Jupiter. While 
it succeeds in producing hot Jupiters that are highly 
inclined with respect to stell ar spin, including on es that 
are retrograde in projection ( Triaud et al.||2010 l, and is 
likely responsible in a number of specific cases (such as 



HD 80 606b N aef et al.||2001| ILaughlin et al.|[20 09; Po^ 
eF^[2009[]Winn et al.||2009| [Hebrard et ai:||20loE| 7 
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it does not preferentially yield low-mass hot Jupitersj^ 
and its effectiveness may be hampered by the pr esence 
of other planets in the system ( Wu fc Murray|2003 ). Fur- 
thermore, population studies establish that only 10% 
of hot Jupiters can be ex plained by Kozai migration 



TABLE 1 

Initial conditions for the example system 



due to binary c ompanions (|Wu et al.||2007| |Fabrycky & 
Tremaine||2007| 

Mechanisms that rely on eccentricity excitation, such 
as Kozai migration or planet-planet scattering, naturally 
produce hot Jupiters that tend to avoid th e region inside 
of twice the Roche radius (Ford & Rasio 2006). However, 
only Kozai migration can naturally explain the 3-day 
pile-up, as the eccentricity rise in this case is gradual and 
planets are accumulated at the right location. The secu- 
lar chaos mechanism described in this paper also leads to 
gradual eccentricity excitation and can therefore inherit 
much of the success of the Kozai mechanism. 



(20081 
itart- 



The noteworthy simulations of Nagasawa et al 
combine planet scattering with Kozai oscillations, 
ing from very compact systems of three equal-mass plan- 
ets, their planets frequently scatter one another onto 
highly inclined orbits, which in some cases triggers Kozai 
oscillations. Their particular set-up yields hot Jupiters 

30% of time, with orbital inclinations that are roughly 
isotropic. The production of hot Jupiters by inter-planet 
Kozai oscillations has also been studied in INaoz et al.l 
(2010). Such a mechanism appears promising if eccen- 



tricities and inclinations can reach order-unity values. 

2.3. Secular Interactions 

Secular interactions are a simplified version of inter- 
planetary interactions, where one can account for the 
forces between two planets by calculating the torque be- 
tween two mass wires. The latter are made by spread- 
ing the mass of a planet along its orbit, weighted by the 
amount of time it spends at that segment. This describes 
the dynamics adequately as long as the planets have no 
close-encounters and do not lie near mean-motion reso- 
nances. Secular interactions allow planets to exchange 
angular momentum but not energy. So planets' semi- 
major axes are unchanged. The long-term evolution of 
the inner Solar system, for instance, is primarily secular 
in nature (Laskar 1989). 

Under certain circumstances, secular interactions can 
gradually raise the eccentricity of the inner planet to near 
unity, even when the initial eccentricity are as low as that 
expected of planets emerging out of dissipative gaseous 
protoplanetary disks. This is the work of chaos. 

When planet eccentricities and inclinations are small, 
secular dynamics is fully described by a linear summa- 
tion of secular eigenmodes that are independently oscil- 
lating, abidin g by the so-called Laplace- Lagrange the- 
ory (see, e.g. Murray & Dermott 20001. The multi- 
periodic variations in the eccentricity of the Earth are 
largely caused by the interference between the eigen- 
modes. These have been claimed to drive climate changes 
( [Milankovitch_1941 ). 

When eccentricities and inclinations rise, linear eigen- 
modes no longer describe the dynamics adequately. Non- 
linear effects can occur. One example is the appear- 

* Kozai migration can readily migrate massive planets inward, 
and perliaps can account for the prese nce of massive hot Jupiters 
found in binary systems | |Zucker &: Mazeli_2002 ). 
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ance of nonlinear secular resonances, including new fixed 
points and separatrices that are not present in the lin- 
ear sy stem (Michtchenko & Malhotra 2004; Michtchenko 
et al. 2006; Migaszewski & Gozdzicwski 2009 K A sec 



ond is the chaotic motio n associated with the overlap of 



neighboring 


resonances (|Sidlichovsky|1990 


Michtchenko 


et al. 112006 
1 \ — " ' 


Lithwick & Wu 2010). 



chaos requires eccentricities and/or inclinations of order 
unity. If the two planets are coplanar, energy and angular 
momentum conservation constrain the motion to be reg- 
ular and quasi-periodic. But if the planets are sufficiently 
inclined, the Kozai resonance can be triggered, leading 
to instability and/or chaos. This requires mutual incli- 
nations > 40° if the initial eccentricities a re very large, 
and ^ 40° for more modest eccentricities (Michtchenko 



et al. 2006 Naoz et al. 20101). For instance, tor the in- 



ner two planets listed in i'able [T] to interact to produce 
ei > 0.98, a mutual inclination of > 85° is necessary. 

For a system of three or more planets, however, the 
threshold for chaos is much reduced. Moreover, the na- 
ture of the secular chaos is different. In the two planet 
case, chaotic systems are quickly unstable on the secular 
timescale, ^10 Myr for typic al parameters ( Michtchenko 
et al.||2006[ |Naoz et al!1|2010[ ). By contrast, in the multi- 
planet case a multitude of resonances can overlap. This 
leads to a gentler type of chaos, with the orbital ele- 
ments gradually diffusing over many secular times, while 
the values of the eccentricities and inclinations remain 
modest. An example of the latter type of behavior is 
the inner Solar system, where chaos is prevalent even at 
eccentricity/inclination levels of a few perce nt, and the 



timescale of the e volution is ^ Gyr ( ,Laskar||1989 Lith- 
w ick fc Wu|[20T0l ). 

In the following, we investigate a planetary system 
with three mildly eccentric and inclined planets to 
demonstrate the appearance of this new type of secu- 
lar chaos — new, that is, in the context of extra-solar 
planetary systems. 

3. SECULAR CHAOS: A WORKED EXAMPLE 
3.1. Numerical Example 

Initial conditions: Parameters for the example sys- 
tem we investigate are listed in Table [T] The inclinations 
are measured relative to the system's invariable plane. 

For the semi-major axes, we space the planets suffi- 
ciently far apart and away from any major mean-motion 
resonances (orbital period ratios are: 1 : 14.7 : 64). The 
choice for the masses is somewhat arbitrary, except for 
our choice that the innermost planet be the least massive, 
which facilitates its excitation. 

The somewhat odd-looking choices for the other orbital 
elements place most of the secular energy (i.e., angular 
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momentum deficit) in the outer planets, or more specif- 
ically in the secular eigenmodes associated with those 
planets (more detail in Q. We find that the occur- 
rence of secular chaos is not particularly sensitive to our 
choices for these values, as long as the system has suffi- 
cient amount of angular momentum deficit (Q. 

Numerics: This simulation was performeausing th e 
SWIFT symplectic integrator ( Levison &: Duncanjl994 ), 
supplemented with routines that model tidal dissipation 
in the planet and in the star, precessions due to gen- 
eral relativistic effects, and precessions due to rotational 
and tidal bulges on both the inner planet and the star. 
Details are presented in the Appendix. These effects are 
essential for determining the final positions of the hot 
Jupiters. 

Numerical precision at extremely high eccentricities is 
of concern. So we adopt a time-step that is 1/100 of 
the inner planet's orbital period (1 yr) for most of the 
integration, but switch it to a value 100 times shorter 
whenever the periastron of the inner planet reaches in- 
ward of 0.1 AU from the star. We find that such a change 
of time-step, even though it breaks the time-symmetry 
of the symplectic integrator, is required to maintain sat- 
isfactory energy and angular momentum conservation. 
The fractional energy error, integrated over an episode of 
extremely high eccentricities (which typically lasts ~ lO"' 
yrs), remains below lO"** as long as ei < 0.98, sufficiently 
small for our problem at hand. 

The angular momentum error is much smaller. 

Results (Figs. Il](2]): Our fiducial 3-planet system 
is chaotic due to secular interactions. The random ex- 
changes of angular momentum (but not energy) among 
planets induce large fiuctuations in their eccentricities 
and inclinations. The orbit of the inner planet, start- 
ing from initial values of ei — 0.07 and ii — 4°, by 
300 Myrs has diffused to ei = 0.985 and h = 70°. It 
has lost almost all of its initial angular momentum to 
the outer planets; equivalently, it has stolen some of the 
outer planets' angular momentum deficit. 

We find that the orbital elements of the inner planet 
undergo a random walk to most of the phase space al- 
lowed by the total energy and angular momentum. How- 
ever, there are a few forbidden regions. The most im- 
portant one is the region of very high eccentricity. The 
inner planet prudently avoids the Roche zone. This oc- 
curs because when the pericentre of the inner planet 
(ai(l — ei)) approaches the star to within a few stel- 
lar radii, finite-size effects (quadrupole precessions asso- 
ciated with the tidal and rotational deformations on the 
planet and the star, respectively), as well as general rel- 
ativity, combine to suppress the secular forcing. Were 
it not for these additional precessional effects, the inner 
planet could be driven to tidal disruption and merger 
with the central star. But as it is, the inner planet stays 
within ei < 0.985, or ai(l — ei) > 0.017 AU in our sys- 
tem. We explain this in the following. 

Due to secular interactions with other planets, the in- 
ner planet's longitude of pericentre precesses at the rate 



is disturbed by the prograde precession induced by the 



(1) 



where Mp is the mass of the perturbing planet, and 
a = ai/ttp is the ratio of the semi-major axes. This 
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1939 


Shakura 


1985 



the following parameters: a Jupiter-like planet with a 
spin period of 3 days, orbiting at a = 1 AU around a 
Sun-like star that is spinning with a period of 10 days, 
we conclude that the tidal bulge raised by the star on 
the planet dominates the precession at high eccentrici- 
ties with the GR effect following not far behind. The 
orbit-averaged precession due to the tidal quadrupole on 
the planet is 



~dt 



I tide 



= 7.5fc2ni 



(1 + 3e?/2 + et/8) 

Ml 



(1-ef) 



2A5 



ai 



(2) 



where k2 is the tidal Love number of the planet, taken 
to be k2 ~ 0.26, and ni, ei, ai and Ri are the planet's 
mean-motion, eccentricity, semi-major axis and radius, 
respectively. Since this rate rises steeply as the planet 
approaches the star, we expect that the secular driving 
is arrested when the planet's pericentre reaches inward 
of 



ai(l-ei)-0.015AC/ 



a 

1/6 



Ml 
Wj 

3/5 



-1/5 



RjJ 



Mp 
Ml 



-1/5 



2/5 



(3) 



This stalling radius is independent of the planet's initial 
position (ai), weakly dependent on the planetary and 
stellar masses, as well as the planet spacing in the system 
(i.e., the value of a). It does, however, scale with the size 
of the planet linearly. Since Jupiter-like planets have a 
fairly uniform radius, the stalling distance spans a narrow 
range for a wide range of system parameters. 

When the orbit of such a planet is tidally circularized, 
the final semi-major axis is moved to twice its stalling 
radius, ~ 2 x 0.015 - 0.03 AUF] This explains why 
the hot Jupiters are piled up at the distance they are 
observed today (Fig. [2]). 

The strength of the dissipative tide can also have an 
affect on the final orbit. Even if the above precessional 
effects are absent, the progression of the inner planet 
toward the star will be halted by the dissipative tide, 
albeit at a somewhat closer distance. 

4. ANALYSIS 

For secularly interacting systems, there is an impor- 
tant conser ved quantity, the angular momentum deficit 
(AMD, e.g. |Laskar|p97 l 



N 
k=l 



AMD - > ;Afc(l- ./l-elcosife), 



(4) 



where N is the number of planets, A^ the circular angular 

momentum of planet k, A^ = ^^^^^ ^G{Mq + mk)ak, 

ik is its inclination relative to the invariable plane (nor- 
mal to the total angular momentum). Since the total 



^ Angular momentum (oc ^ a(l — e'^)) is roughly conserved dur- 
ing tidal dissipation. So the post-circularized a'j^ is related to the 
pre-circularized ai via a'j = ai(l — e^) ~ 2ai(l — ei). 
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Fig. 1. — Formation of a hot Jupiter in our fiducial system (SWIFT integration with tides and GR). Left: radial excursions of the three 
planets (semi-major axis, periapse and apoapse) are shown as functions of time, with the various radii relevant for hot Jupiters marked by 
arrows; right: planet inclinations measured relative to the system's invariable plane. All planets initially have mildly eccentric and inclined 
orbits, but over a period of 300 Myrs, so much of the angular momentum in the inner-most planet can be removed that its eccentricity and 
inclination can diffusively reach order unity values. Planet interactions leave all semi-major axes unchanged, a tell-tale sign that secular 
interactions dominate the dynamics. At ~ 300 Myrs, the pericentre of the inner planet reaches inward of a few stellar radii and tidal 
interaction with the central star kicks in (details in Fig. [2jl. Precessions by general relativity, by tidal and rotational quadruples, as well 
as tidal dissipation, prevent the pericentre from reaching mward of the Roche radius. As a result, the final hot Jupiter has a period of ~ 3 
days. 




2.95x108 3xl0« 2.95xl0« 3xl0« 

time (yrs) time (yrs) 

Fig. 2. — Same as that in Fig. [l]but expanding the time axis around 300 Myrs to highlight the process of tidal circularization. Secular 
chaos raises the inner planet's eccentricity diffusively to a maximum value of 0.985, and decreases its periapse to ai * (1 — ei) ^ 0.015 AU 
- as determined by a balance between secular forcing and close-range forces. When this occurs (at time 294 Myrs), tidal dissipation in the 
inner planet removes orbital energy while conserving orbital angular momentum. This brings the planet from an initial orbit of ai = 1 AU 
to an orbit of ai (1 — 0.985'^) X 1 AU ~ 0.027 AU. It is then dynamically decoupled from the outer two planets. Since the total angular 
momentum deficit (AMD) in the system is absorbed by the inner planet (and subsequently removed by tidal dissipation), the outer two 
planets lose AMD, and hence become more circular and coplanar after the hot Jupiter has formed. The inclination of the final inner orbit 
relative to the invariable plane is ~ 70° (right panel). The evolution is highly chaotic, and a slight modification of the initial conditions 
changes the evolution dramatically (see Fig. [sl. 
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angular momentum is conserved, and secular interac- 
tions do not modify the orbital energies {ak constant), 
the AMD is conserved during secular interactions. For 
circular, coplanar systems, the AMD is zero, and the 
AMD increases with increasing e's and Vs. 

Here, we analyze our numerical results to illustrate the 
condition for hot Jupiter formation. We find that a suf- 
ficient amount of AMD is requisite. Firstly, the value 
of AMD limits the maximum eccentri city and inclina- 
tion an individual planet can attain ([4.1 1. Secondly, 
only when AMD is large enough, can it be snared among 
different secular eigenmodes (equipartition), ultimately 
driving the inner planet to extreme orbits. We call this 
sharing process 'secular chaos'. We illustrate the deep 
analogy between AM D an d kinetic energy in a thermo- 
dynaniical system in §4.2[ and analyze the diffusi ve an d 
cha otic nature of the energy sharing process in §4.3| & 
|4.4[ We also briefly look at the issue of AM D ge neration 
by mean-motion resonances in the system (^4.5). 



4.1. Maximum Eccentricity and Inclination 

To be propelled into the hot Jupiter status, the inner 
planet has to reach so close to the star that tidal dissi- 
pation operates. Let this be roughly ai(l — ei) < 0.05 
AU. If all AMD can be transferred to the inner planet, 
this condition translates into 



AMD > Ai 



1 



ai 



O.IAU 



-1/2 



COS«i 



(5) 



So a planet that is closer to the star and lower in mass 
needs less AMD to become a hot Jupiter. From now on 
we measure AMD in units of the circular angular momen- 
tum of the inner planet (Ai — 1). The above condition 
translates into 



AMD > 1 - 0.3 cos ii 



(6) 



for fli = 1 AU. So to produce a coplanar hot Jupiter (ii = 
0), AMD > 0.7, while a retrograde hot Jupiter would 
require AMD > 1.0. Fortunately, this is not difficult to 
satisfy - the outer planets can carry plenty of AMD even 
at low values of eccentricity/inclination. Our example 
system has AMD = 1.17. Retrograde hot Jupiters can 
potentially be produced. 

4.2. AMD and Kinetic Energy 

The AMD is for a secularly interacting system what 
kinetic energy (or temperature) is for a thermodynamical 
system. This analogy runs deep, as we shall show here. 

We intro duce the coniplex Poincare variables Z k and 
Cfc (see, e.g. |Laskar||1997| [Murray fc Dermott]|2000| ), 



Zk -■ 



V2\ 1 - a/i - elexp{in7k), 



Cfc = V2W W 1 - 6^(1 - coszfc) exp(ir2fc). 



(7) 



where Wk is the longitude of periapse, and flk the lon- 
gitude of the ascending node. At low eccentricities and 
inclinations, Zk ~ ekexp{izuk), a nd Ck y ik^ xpiiftk)- 
The AMD may then be recast as ( |Laskar||1997[ ) 



N 



Ak 



.4Mi? = ^^(|zfc|^ + |Cfc 



(8) 



The resemblance of AMD to kinetic energy becomes ob- 
vious in this form: while the "inertial mass" for each 
planet corresponds to its circular angular momentum A^, 
the (2:fc,Cfe) P&ir corresponds to its "velocity." 

When AMD is zero, the system will remain coplanar 
and circular and stable forever. At low AMD, secular 
interactions lead to periodic exchanges of angular mo- 
mentum between planets. This, however, is not related 
to the equipartition process. The dynamics can be de- 
composed into that of linear secular eigenmodes. And 
the periodic variations in orbital elements are caused 
by the interference between these modes (the so-called 
Laplace-Lagrange theory) . Each linear mode oscillates at 
its characteristic eigenfrequency with constant amplitude 
and phase. Let the eigenvectors be Zka and Cka- They 
are orthonormal after each component is pre-multiplied 

by v/Afc/2, i.e., J^k ^ZkaZkp = J2k = ^ap, 

where 5aj3 is the Kronecker delta; the overall phase of 
each eigenvector is chosen so that all components are 
real. Projecting the complex orbital elements onto these 

eigenvectors, Zk = Yfl=i^<^Zka,C,k = Lq^at+i ""Cfca^ 
we can re-express AMD as 



2N-1 



AMD = ^ 



(9) 



a=l 



Each eigenmode resembles one degree of freedom in a 
thermodynamical system, and the AMD resembles the 
total energy. In the linear solution, if one mode is initially 
assigned all the AMD, it will retain it forever. As a 
result, each planetary orbit moves within a certain bound 
as given by the initial condition. 

As AMD rises, energy transfer (or really, AMD trans- 
fer) becomes non-periodic and chaotic. Orbital elements 
are allowed to wander as in a random-walk diffusion. 
This may ultimately lead to AMD equipartition between 
different secular modes (different degrees of freedom in 
the system), as well as approximate AMD equipartition 
between different planets ^ The least massive or the clos- 
est planet has the smallest inertia. AMD equipartition 
implies that such a planet can reach very high eccentric- 
ity and/or inclination, providing the condition for hot 
Jupiter formation. 

4.3. Diffusion of AMD observed 

Diffusion of energy in a weakly nonlinear system is 
an extensively studied subject, start ing from the famou s 
Fermi-Pasta-Ulam (FPU) problem (Fermi et al. 1955). 
In the following, we present evidence tor AMD diitusion 
in our example system, and illustrate the criterion for 
AMD diffusion. 

One line of evidence comes from the amplitudes of the 
secular eigenmodes]^ The initial conditions chosen in 
Table [T] correspond to deliberately depositing almost all 

^ There is a trivial inclination mode with zero frequency, which 
corresponds to the overall tilt of the reference plane. Here, we take 
the reference plane to be the invariable plane, so the amplitude of 
this mode is zero. 

^ Equipartition of AMD amongst planets is only approximately 
correct. As an example, the four terrestrial planets in the Solar 
system have similar AMD today. They have likely undergone ex- 
tensive chaotic diffusion in the past. 

* Here, we decompose using the linear eigenvectors even though 
they are invalid at large amplitudes. A more rigorous approach may 
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TABLE 2 

Frequencies and Amplitudes (initial and final) of the five 

LINEAR secular EIGENMODES. 



mode 


frequency 


amplitude (|aa|) 




(arcsec/yr) 


t = 


2.9 X 10* yrs 


el 


4.83 


0.03 


0.63 


e2 


8.30 


0.31 


0.21 


e3 


2.04 


0.83 


0.51 


il 


-4.61 


0.14 


0.24 


i2 


-10.57 


0.62 


0.64 



AMD into the secular modes associated with the outer 
planets, and little in the eccentricity mode associated 
with the inner planet. There are three advantages to 
this. First, AMD transfer between modes occurs on sec- 
ular or longer timescales. So for the first tens of millions 
of years, all three planets have small e's and i's, as appro- 
priate for planets emerging from a protoplanetary disk. 
Second, by initializing the inner eccentricity mode with 
low amplitude, we are furthest away from our preferred 
end state, when this mode acquires enough AMD to place 
the inner planet on a e ~ 1 orbit and a tidal encounter 
with the star. Third, by concentrating AMD into only 
a few modes, we can best observe the approach towards 
equipartition. 

Table [2] and Figure |3] show that after 294 Myrs of evo- 
lution (right before tidal dissipation sets in) , the AMD of 
mode 1 has grown diffusively by a factor of ^ 20^ ^ 400, 
while the total AMD — J^a Wa]^ remains constant. The 
normalization of \aa\ in Table ^ is such that the circular 
angular momentum for the inner most planet (at 1 AU) 
has the numerical value of one. AMD equipartition is 
reached in ^ 2.5 x 10^ yrs. 

A second way to quantify diffusion is by tracing the 
"spectral entropy" ( |Livi et al.||1985| [Goedde et al.||1992 | 



S 



2N-1 

Ea In Ea , 



(10) 



where Ea = laap/X^akaP is the fractional energy 
(AMD) in mode a. This entropy increases from its min- 
imum value of (when one mode has all the AMD) to 
the maximum value of InA^mod = lii(2iVpi — 1) = ln5 as 
the system difi'uses toward energy equipartition. Spectral 
entropy is analogous to entropy in a thermodynamical 
system. 

Fig. [3] shows that for our example system, AMD is 
gradually shared among different eccentricity and incli- 
nation modes, and the spectral entropy rises toward In 5 
over a 10* year timescale. 

A third way to observe diffusion is to plot the probabil- 
ity distribution function of ei , following Laskar ( 2008 ) , as 
shown in Fig. [4] For our example system, ei is broadly 
distributed from to 1, wit h its distribu tion roughly de- 
scribed by a Rice function ( Laskar|[2008 ) , 



■ exp 



2(72 



(11) 



where Iq is the modified Bessel function of the first kind 



be to decompose using nonlinear eigenvectors, as is done in ( |Laskar| 
|2008||. But our approach suffices for the purpose of illustrating 
TCSTD diffusion. 



with order 0. This is the expected distribution for \z\ if 
z = x + iy with x and y being two independent Gaussian 
variables with mean v and variance a. 

What is required for AMD diffusion to occur? We find 
that there must be a sufficient amount of AMD in the 
system. For instance, when integrating the example sys- 
tem, but with all secular mode amplitudes reduced by 
a factor of 2 (total AMD reduced to 1.17/4 = 0.29, the 
'Twice-Lower' curve in Fig. El), or when flattening all 
orbits into coplanar ones (AMD= 0.72, the 'Coplanar' 
curve in Fig. 4]), we find that AMD diffusion is largely 
suppressed. Wnile AMD allows ei to reach 0.58 ('Twice- 
Lower') and 0.955 ('Coplanar'), diffusion is only able to 
bring ei to 0.15 and 0.58, respectively. Motion is either 
largely quasi-periodic or weakly chaotic. This is in con- 
trast to the example case where the inner planet explores 
phase space and reaches its maximum eccentricities and 
inclinations. 

Why does the amount of AMD make a difference? The 
interested readers are referred to Lithwick & Wu (20101 
for detailed quantitative analysis. Here we only comment 
that to allow AMD diffusion, chaos is essential, and chaos 
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onances. Their widths increase sharply with mode am- 
plitude. A lowering of the mode amplitudes by a mere 
factor of 2 can shrink the resonance width by a large fac- 
tor. This qualitatively explains why the twice-lower case 
is not chaotic. In a coplanar system, all resonances that 
involve the inclination modes are ineffective, and with 
fewer resonances to drive chaos, they dynamics become 
much more regular. 

4.4. Chaotic Processes Leading to Hot Jupiters 

Given the chaotic nature of the dynamics, the final 
state of the system depends sensitively on initial con- 
ditions. In our example case, we obtain a hot Jupiter 
with ai = 0.027 AU and an inclination of 70° at the 
end of the 300 Myrs integration (Fig. [2]). In another 
case, modified from the previous one in energy and an- 
gular momentum by about 0.1% at time 293 Myrs (just 
before the onset of tidal dissipation), we find that the 
inner planet narrowly avoids being tidally circularized 
into a hot Jupiter straight-away, and is able to return 
its AMD to the outer planets. The eventual hot Jupiter 
thus formed has ai = 0.04 AU and an inclination of 25° 
(Fig. pi). Our other experiments show that in some 
cases tne tidal dissipation process is gradual and occurs 
over many episodes of high eccentricities. There could 
be an extended period during which the planet destined 
to become a hot Jupiter is temporarily parked at an in- 
termediate distance (e.g., oi =0.1 AU) with large eccen- 
tricities and inclinations. They may help to explain the 
presence of 'warm Jupiters'. 

With a tidal quality factor for the planet of Qp = 10^, 
the example case sees most of the orbital binding energy 
deposited inside the planet within a couple Myrs, with 
an averaged (over 1 Myrs timescale) heating rate of 5 x 
lO^^ergs/s, ~ 15 times higher than the self-luminosity 
of a Teff = 1000 K Jovian planet. The total tidal energy 
deposited is 1.5 x 10^^ ergs, again about 15 times larger 
than the gravitational binding energy of the planet. This 
has the potential of disrupting the planet unless the heat 
is deposited in regions of short thermal time. However, if 



8 




IQS 2x109 3x10^ 108 2x108 SxlO^ 

t (yrs) t (yrs) 

Fig. 3. — The diffusive behaviour for the system in Fig. IT] demonstrated using eigenmode energies (really, AMD's) and spectral entropy, 
The curves in the left panel are, from top to bottom, total AMD, AMD in eccentricity modes, AMD in inclination modes, and AMD in 
the eccentricity mode associated with the inner planet, all plotted as their ratio to the total AMD (which is conserved to better than a 
percent before tidal dissipation sets in). For clarity, the last three are plotted as running overages (over 2 X 10* yrs). The first eccentricity 
mode gains AMD from other modes, and there is AMD exchange b etwe en eccentricity and inclination modes. The right hand panel shows 
the exponential of the spectral entropy, exp(S'), as expressed in eq. |10[ The process of AMD equipartition raises exp(S') from 2 (2 modes 
dominating the AMD ) to ~ 5 (all modes share comparable AMD, marked by the dotted curve). When tidal dissipation sets in, the total 
AMD decreases, leading to a drop in the entropy. If we turn off the dissipation and integrate the system further, we observe that the 
spectral entropy fluctuates, but it seldom departs from the equipartition value for an extended period. 
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Fig. 4. — The probability distribution of the inner planet's eccen- 
tricity, dN/dei, in the 'Example' system of Fig. [l] in a system with 
twice lower eccentricities and inclinations ('Twice-Lower'), and in 
a system that has the same eccentricities but is coplanar ('Copla- 
nar'). The first case is integrated for 300 Myrs (stopped short when 
the planet reaches the zone of tidal circularization) , while the latter 
two are integrated for 1.5 Gyr. The maximum eccentricity allow- 
able in each case (see Fig. [7| are marked by the arrows. Diffusion 
in the latter two cases is partially inhibited. The dashed line is a 
Rice distribution with mean 0.2 and variance 0.25. The example 
system is expected to approach the Rice distribution after a long 
time. 

the tidal dissipation process is more gradual, the mode of 
heat deposition can be drastically different, potentially 
impacting on the final sizes of hot Jupiters. 

4.5. Effects of Mean-motion Resonances 

Mean- motion resonances (MMRs), not to be confused 
with secular resonances, can also affect secular chaos. 



This is true even if the MMR's are of high order. Because 
MMRs can change the semi-major axes, the dynamical 
system has more degrees of freedom and can explore dif- 
ferent parts of phase space, potentially enhancing diffu- 
sion. We experiment by placing the outer two planets 
near the 2 : 1 MMR (6 AU and 9.52AU). We find that 
diffusion proceeds quickly even with an AMD as low as 
0.29, a value that corresponds to our twice- lower system 
(which shows regular behaviour). But when we move 
the outer planet outward by a mere 0.5 AU, the system 
behaves regularly again. This confirms that, at least in 
this case, MMR's can be responsible for facilitating AMD 
equipartition. 

In addition, even for systems which are not initially 
next to any lower order MMR's (such as our example 
case), as eccentricity and inclination rise to order unity 
values, low order MMR's are activated. For the case 
shown in Fig. [ij the semi-major axis of the inner planet 
undergoes increasingly large variations at late stages of 
the evolution, perhaps indicating MMR's at work. Ac- 
companying this is the non-conservation of AMD. Are 
MMR's a significant source of AMD? If so, they could 
allow the inner planet to reach a higher eccentricity than 
allowed by the initial AMD of the system. This intrigu- 
ing possibility deserves to be explored. 

Another possible connection between secular chaos and 
MMR's exists. In systems which are initially more com- 
pact than our example system, but not compact enough 
to have immediate close encounters, secular chaos in- 
creases the eccentricity of inner planets, allowing MMR's 
to function at later stages, leading finally to planet scat- 
tering well after the protoplanetary disks have dissipated. 

5. PREDICTIONS FOR HOT JUPITERS 

We have demonstrated that secular chaos can produce 
hot Jupiters. The remaining central issue is how preva- 
lent this mechanism is. Is it prevalent enough to explain 
the observed frequency of hot Jupiters? 
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Fig. 5. — Evolution of a system that is slightly modified from that in Fig. [2]by 0.1% in total orbital energy and angular momentum, at 
time t = 2.93 X 10* yrs. The subsequent trajectory differs. The inner planet in this case suffers some tidal dissipation at t = 294 Myrs, 
and is migrated to ai = 0.75 AU. However, it avoids being turned into a hot Jupiter right-away, and returns its AMD to the outer planets 
(left panel). Secular chaos continues to operate until the inner planet is turned into a hot Jupiter with ai = 0.04 AU and an inclination of 
25° . This demonstrates the sensitive dependence on initial conditions. 



Unfortunately, since we do not know the initial config- 
urations of planetary systems, it is difficult to predict the 
frequency of hot Jupiter production by secular chaos. We 
also face the problem that a systematic survey of the rel- 
evant parameter spaces is numerically expensive, at least 
using our current technique of N-body integration. So in 
the following, we discuss qualitative predictions based on 
our present understanding of secular chaos. 

5.1. General Predictions 

The predictions are ranked roughly in decreasing order 
of certainty. 

1. a pile-up of hot Jupiters around three-day orbital 
periods. 

The characteristic stalling at 3-day orbital periods 
that we observe in our simulations is explained by 
a combination of tidal precession and tidal dissi- 
pation (both due to tides on the planet). Preces- 
sions by other close-range forces (GR, rotational 
quadrupole) are less important. The tidal preces- 
sion stalls the rise of the eccentricity to about a few 
times the Roche radii, and tidal dissipation finishes 
the job by circularizing the highly eccentric orbit 
to that of a hot Jupiter. Location of the pile-up 
remains largely unchanged even if the rate of tidal 
dissipation is orders of magnitude weaker. This 
conclusion also applies to Kozai migration, another 
secular migration process. 

2. hot Jupiters are lower in mass compared to other 
giant planets. 

This prediction stems from AMD conservation. To 
become a hot Jupiter, the inner planet has to reach 
an eccentricity so high that ai(l — ei) < 0.05 AU. 
For a given amount of AMD, a lower mass inner 



planet can reach a higher eccentricity. Compared 
to giant planets at larger distanc es, there is a clear 
deficit of massive hot Jupi ters (Zucker & Mazeh 
[2002; ,Udry fc Santos|[2007| ). Ko zai migration, in 
contrast, has no mass preferences ( Wu et al.|2007 ). 



3. hot Jupiters have no companion within a few AU, 
but have companions roaming at larger distances. 

On the one hand, the high eccentricity episode 
the inner planet undergoes before it is tidally cap- 
tured implies that hot Jupiters have no compan- 
ions within a few AU. This agrees with observa- 
tions where hot Jupiters appear to be strikingly 
alone ( Wright et al. 2009 1 , and where attempts at 
measuring transit timing variations have repeat- 
edly turned up empty-handed. On the other hand, 
secular chaos requires driving by other giant plan- 
ets. They should be roaming at large distances 
(outside a few AU) and remain to be detected - 
some may have already shown up a s radial veloc- 
ity residuals in hot Jupite r systems (Fischer et al 



2001 Wright et al. 



2009). There is also the in- 



triguing possibility tnat one of the companions is a 
binary star. 

4. frequency of hot Jupiters should rise with stellar 
age. 

Being a diffusive process, secular chaos operates on 
timescales comparable to or longer than the secular 
precession timescale. The latter, for our fiducial 
system, is of order M^/M2Pi/Pi ~ 10^ yrs. As 
such, we expect stars that are within a few tens of 
million years after their disk dispersal should have 
a lower hot Jupiter fraction than stars that are a 
few Gyrs old. This prediction should be quantified 
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by extensive numerical simulations that start with 
reasonably realistic planetary configurations. 

5. orbits of hot Jupiters could be strongly misaligned 
with stellar spin. 

We will discuss the planet inclination in j |5.2| 

If secular chaos is responsible for producing hot 
Jupiters, this has a number of implications for planetary 
systems at large. 

1. a large number of planetary systems should have 3 
or more giant planets on (mildly) eccentric, inclined 
orbits. 

Only systems that have sufficient AMD can make 
hot Jupiters. Hot Jupiters may be the tip of the 
iceberg in terms of their system AMD values. If so, 
then most giant planets we know of should reside 
in systems with three or more giant planets. There 
should still be residual AMD in the outer planets. 

2. warm Jupiters could arise from secular chaos and 
their orbit could be misaligned with respect to the 
stellar spin. 

We call giant planets at a few tenths of an AU 
'warm Jupiters,' for the reason that they are at in- 
termediate locations between hot Jupiters (~ 10~^ 
AU) and cold Jupiters a couple AU). 

Among those that fail to make hot Jupiters, there 
may be warm Jupiters - planets which have suf- 
fered some degree of tidal dissipation but have yet 
to be tidally captured into hot Jupiters. These 
planets are temporarily intercepted from their in- 
ward spiral by interactions with outer planets that 
can inject into their orbits a fresh boost of angu- 
lar momentum. Eventually, these planets would be 
dragged in to become hot Jupiters, but while they 
are in their temporary parking space, their orbits 
could have high inclinations as well as high eccen- 
tricities. Observationally, there is a 'period val- 
ley' at these distances, indicating that perhaps the 
warm Jupiter phase is relatively short-lived. Mea- 
surements of spin-orbit angle for these planets will 
be useful constraints. 

3. even for systems where the inner planets can not 
reach high enough eccentricity to be captured into 
hot Jupiters, secular chaos may have observable 
consequences. 

Starting from mildly eccentric, inclined orbits, in- 
ner planets in planetary systems may gradually 
extract AMD from outer planets. The long-term 
eccentricity and inclination distributions for these 
planets approach the Rice distribution (Fig. [4]). 
This reduces to a Rayleigh distribution when the 
centroid is much smaller than the variance. Fig. [6] 
shows the observed eccentricity distribution of Jo- 
vian planets: it may be represented by a Rayleigh 
distribution with variance 0.25. If future Rossiter- 
McLaughlin measurements provide us with a simi- 
lar distribution for the orbital inclinations, this will 
be a strong proof in favour of secular chaos. 
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Fig, 6, — The eccentricity distribution of Jovian planets that 
have pericentre between 0.1 and 10 AU and projected masses above 
O.lMj. (the solid histogram). The lower limit of 0.1 AU is selected 
to preclude planets that have undergone significant tidal circular- 
ization. The dotted line is the Rayleigh distribution with a = 0.25. 
The sub-population of more massive planets with M s'mi > 4Mj 
are represented by the dashed histogram. It appears that more 
massive planets have a hotter distribution. 



If secular chaos is largely responsible for the ob- 
served planet eccentricities, one expects that, for 
the same amount of AMD available in the system, 
lighter planets in general reach higher eccentrici- 
ties. However, data (Fig. p| portray an inverse 
correlation - more massive [Alsini > AMj) plan- 
ets have a hotter eccentricity distribution. This 
may be interpreted as due to pollution from sys- 
tems undergoing Kozai cycles induced by stellar 
companions, which is capable of exciting eccentric- 
ities even in very massive planets. This interpreta- 
tion requires that massive planets occur preferen- 
tially in binary systems. 

4. secular chaos can stabilize planetary systems. 

In our simulation (Fig. IT]), the tidally captured hot 
Jupiter removes AMD from the system, stabiliz- 
ing the outer planetary systems as a result. This 
could be a generic process in organizing planetary 
systems on long timescales. 

5. can secular chaos explain hot Neptunes or hot 
Earths? 

If the inner planet has a lower mass, it experi- 
ences a weaker near-range precession due to its 
smaller tidal and J2 moments. Eq. |3] predicts 
that a Neptune-like planet (M = M^r = l/17Mj, 
and R = Rn = 0.36i?j), driven to chaos by secu- 
lar forcing of giant planets at a few AU, could be 
stalled at periapse distance ~ 0.009^[/, and hence 
it would be circularized at ~ 0.018AC/ ~ 4i?0. 
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The final inclination of the inner planet deserves an 
i n-depth discussion 



Fig. 7. — The eccentricity-inclination phase space that is tra- 
versed by the same system as Fig. [ijover a total of 400 Myr, in 
the absence of tidal dissipation. The thick red curve bounds the 
region within which the inner planet can explore, as dictated by 
energy and angular momentum conservation (assuming constant 
ai, so this is determined by the initial AMD). Secular chaos trans- 
ports the inner planet to very high eccentricities (bounded at the 
right by near-range precessions), and to retrograde orbits. If tidal 
dissipation were turned on, it would operate rightward of the ver- 
tical line at e = 0.95 (pericentre distance of ai(l — ei) = 0.05 
AU). The right panel shows the distribution of inclinations when- 
ever ei > 0.95. This can be regarded as an approximation to the 
distribution of final inclinations. For this system, retrograde orbits 
are achieved a few percent of the time. 



This distance would be approximately doubled if 
it is forced by other Neptune-mass planets, while 
a more inflated planet would also be stalled at a 
greater distance. Hot Neptunes thus produced will 
tend to have a broader pile-up than hot Jupiters, 
because the masses of the perturbers can extend 
over a large range. Such a mechanism could ex- 
plain some of the observed hot Neptunes - for in- 
stance, the HD 125612 system where a hot Neptune 
is accompanied b y two Jupiter-like plan ets at much 
larger distances ( [Lo Curto et al.pOlO ). 

Bouchy et al. ( 2009 1 noted that ~ 70% of low mass 
close-m planets (with Msini < O.lMj) have de- 
tected planetary companions. This differs from the 



hot Jupiter case. It likely relates to the lower re- 
quirement for the planetary spacing when the inner 
planet is less massive. 

An Earth-like planet, on the other hand, will not be 
stalled at large enough distances to avoid Roche- 
lobe overflow. Secular chaos is an unlikely agent 
for making hot earths, especially considering that 
the tidal damping time for an Earth-like planet at 
0.03 AU exceeds a Hubble time. 

5.2. Stellar Obliquity 



Fabrycky & Winn ( 2009 1 found that a majority of tran- 
siting planets have their orbits aligned with the stellar 
spin. However, that majority was quickly weak ened by 



a slew of misaligned systems (Triaud et al. 2010j ). The 



current situation is unclear, partly due to small number 
statistics, though a correlation between stel lar obliquity 
and s tellar spectral type has been noted (Winn et al. 
20T0| ). 

Cmr example system (Fig. |1|) produces a hot Jupiter 



with an orbit inclined by 70° from the invariable plane of 
the planetary system, while its near-identical twin yields 
a hot Jupiter with an inclination of 25° (Fig. [5| . Secular 
chaos could in principle drive the inner planet to inclina- 
tions between and 180° relative to the invariable plane 
(Fig. [T]), given the amount of AMD in this system. To 
demonstrate this, we turn-off tidal dissipation in our ex- 
ample integration, and evolve the fiducial system for 400 
Myr. The eccentricity and inclination phase-space tra- 
versed by the inner planet is presented in Fig. [Tj Within 
the constraint of AMD conservation, the planet is able 
to diffusively reach high eccentricity, with both prograde 
and retrograde orbits. The final inclination for the hot 
Jupiter then depends on chance. 

In Fig. [TJ we also display the distribution of inclina- 
tions whenever 1 — ei < 0.05, as a proxy for the final 
inclinations of hot Jupiters. We find that retrograde or- 
bits represent a few percent of the total. The preponder- 
ance of prograde orbits is related to the prograde initial 
condition. 

This result is specific to our example system and is 
relative to the invariable plane. The invariable plane 
does not necessarily coincide with the stellar equatorial 
plane ~ the relative angle is only 7° in the Solar system, 
but it can be large elsewhere]^ 

6. SUMMARY 

Hot Jupiters, while representing only a small fraction 
of all known extra-solar planets, demand special atten- 
tion. They are most at odds with planet formation the- 
ory; they are detected disproportionately in radial ve- 
locity and transit surveys; and they are most accessi- 
ble to characterization. Their rarity may indicate that 
their formation requires extreme circumstances. How- 
ever, they may teach us much about the general condi- 
tions of planetary systems. 

Hot Jupiters are piled up around 3-day orbital peri- 
ods, with rapid cut-offs both inward and outward of this 
distance; they tend to be less massive than more distant 
planets; many of them have orbits that are misaligned 
relative to the stellar spin; and they are remarkably anti- 
social in having few detected companions. 

In this work, we show that most of these characteristics 
can be explained if hot Jupiters are produced by secu- 
lar chaos. These planets, originally located at > 1 AU, 
acquire AMD from planets that are further out in the sys- 
tem. The outer planets can be mildly eccentric and/or 
inclined. But the same AMD produces much greater ec- 
centricities and inclinations when it is transported to an 

^ Here, we have assumed that the stellar spin is aligned with the 
invariable normal when calculating precession due to the stellar 
rotational bulge. 
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inner planet, especially if the inner planet is less mas- 
sive than the outer ones. The extreme high eccentricity 
allows the inner planet to reach inward of a few stellar 
radii and be tidally ensnared by the central star into a 
hot Jupiterj^ We find that the criterion for hot Jupiter 
production is a sufficient amount of AMD. 

Only 5 hot Jupiters have known planetary compan- 
ions, and these companions are typically eccentric and 
may have contributed to secular chaos. We note that the 
hot Jupiter (at 0.06 AU) in the Ups And system, where 
two other massive planets orbit at 0.8 and 2.5 AU, with 
eccentricities 0.2 and 0.3 respectively, may well be pro- 
duced by the secular chaos presented here. This possi- 
bility is further boosted by the recent finding that the 
outer two planets ' inclinations are misa ligned by ~ 30° 
from each other (McArthur et al. 2010 ). Similarly, the 
retrograde hot Jupiter m WAyF-8b fQueloz et al.l 2010 1 



is in a stellar system including an M-star companion at 
^ 600 AU. Moreover, the radial velocity trend indicates 
a companion at distance > lAU that is more massive 
than 2Mj. This could also be a hot Jupiter produced by 
secular chaos, with the M-star acting as the third planet. 

Hot Neptunes may also be formed via secular migra- 
tion. However, hot Earths are different . It is not that 
Earth-like planets could not undergo secular chaos, but 
that once they do, they cannot be stalled at a safe enough 
distance to avoid being swallowed by the central star. 

Secular chaos has also been found to be responsible for 



instability in the inner Solar System ( |Laskar||2008[ |Lith , 
wick fc Wu||2010 ). We speculate that secular chaos may 
be a frequent phenomenon in planetary systems. It may 
help to excite inner planets to higher eccentricities or in- 
clinations. If these planets are removed, the remaining 
planetary systems may be stabilized for a time compara- 
ble to the system age. 

The success of this theory depends on two unknown 
factors. One is the amount of initial AMD in the sys- 
tem. The other is the typical configuration of planetary 
systems when emerging out of the protoplanetary disk. It 
also needs to be demonstrated solidly that secular chaos 
can lead to a large fraction of retrograde hot Jupiters. 
Observationally, if not only hot planets, but also warm 
or cold planets can be shown to have significant orbital 
inclinations relative to the spin of their host stars, this 
would boost the case for the ubiquity of secular chaos. 
Future Rossiter-McLaughlin measurements should be ex- 
tended to transiting planets at large distances. 
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performed. YW further acknowledges the NSERC fund- 
ing and a Caltech sabbatical stay. Lastly, we thank M. 
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APPENDIX 

The standard SWIFT package distribution does not contain a treatment for the general relativistic precession, or 
precession due to tidal and rotational bulges. These dominate over secular precession by other planets when the inner 
planet reaches very close to the host star. Together with tidal dissipation, these effects determine the final orbit of 
the planet, as well as the timescale for tidal circularization. Below is our implementation of all these processes in the 
SWIFT code. Since non-Keplerian accelerations in the symplectic SWIFT integrator are incorporated as velocity kicks 
between Keplerian drifts, we need expressions for the perturbative accelerations. 

To first order in (w/c)^, the GR effect can be written as a perturbation to the Newtonian gravitational potential as 

$GR = — ^^T' (1) 

where r is the radial distance between the star and the planet, L = mp^yGM^a{l — e^) is the orbital angular momen- 
tum. The purely radial force associated with this potential, 

„^ 3G^M^a(l - e2) , 
acR = - V<I>GR = f , (2) 

gives rise to a precession of the eccentricity vector. We confir m that such a numerical procedure yields the following 
orbit-averaged precession rate for the longitude of pericentre ( Einsteiii]|1916 | 
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For th e tidal effects, we first consider the rotational bulge and the tidal bulge on the star. The following expressions 
are from Sterne (1939), differing only in notation. Let the stellar spin rate be w*, radius i?*. The centrifugal potential 
due to the stellar spin, and the quadrupole tidal potential due to the planet acting on a point at distance D away from 
the stellar center is 



acting = +^^.'i?'P2(cOS0')-^ 

3 r 



(4) 



Here, the Legendre function P2{x) = 1/2 (3a;^ — 1), and the connecting angles are defined as cos6' = D-f, cos^' = D-cJ*, 
where r is the vector connecting the two bodies. The global distortion of the star under the above potential casts a 
response potential in its surrounding (measured at position D) 



L>3 



^P2(COS0')-^^2(COS( 

3 r'^ 



where ^2 is the Love number and is taken to be 0.029 for the star (polytrope n 
n = l)j^ The barycentric acceleration the reduced particle feels is therefore 



a = 



response; 



(5) 



3) and 0.52 for the planet (polytrope 



(6) 



where the reduced mass /i = m*mp/(m* + nip). 

We trace the planetary motion in the frame of the invariable plane. And we assume that the stellar spin coincides with 
the normal of the invariable plane. So at the heliocentric position of the planet (x, y, z), cos = 1, and cos 6' = jr^ , 



where r = 



Here the z term in the potential gives rise to spin-orbit coupling and the orbit normal of 



the planet precesses around the spin direction 

The bulges on the planet are treated similarly, except we assume that the planet spin is aligned with the orbit 
normal, so cos & = 0. 

The above close-range accelerations are expressed for the barycentric movement. However, one should correct for 
the fact that the kick asked by the SWIFT code is the heliocentric value. 
The dissipative part of the tidal effect is handled in a way that differs from the standard treatment of 



Lee & Peale 



(2003). We use the weak friction prescription for the equilibrium tidal bulge, and calculate the effect of dissipation in 
both the star and the planet. The tidal bulge raised on either body produces a potential given by the second term 
in eq. [5] Due to finite dissipation inside the body, there is a delay between the response and the forcing. We can 
assume either a constant time lag (r), or a constant phase lag (e). In the latter case, we can introduce a tidal Q factor 
(Goldreich & Soter|l966), which is related to the phase lag and the time lag as e = l/Q, and t = ^ 1^^^ respectively. 
Here wtido is the tidal forcing frequency. For the eccentricity tide, wtidc is simply 2n. The acceleration associated with 
the delayed tidal bulge raised on body M (with radius R) by body m at distance r is ( )Hut|,1981^ 



Gm2 fR 



r 
9-r 

r 



f - (w - e)Te 



(7) 



where w is the rotational velocity and 9 the instantaneous orbital angular velocity. In particular, f = nae sin //\/l — 
for a Keplerian ellipse. The first half of the radial term contributes to orbital precession (dealt with above) but no 
energy dissipation (as it is anti-symmetric within a Keplerian ellipse and cancels out over an orbit). We ignore the 
angular force which transfers angular momentum - the spin angular momentum of the planet is much smaller than its 
orbital angular momentum and so we assume that the planet is quickly synchronized with the orbit, while we assume 
that the tidal Q factor associated with the star is so large that there is no angular momentum being transported 
between the orbit and the stellar spin. So we are left only with the second half of the radial force. This is easily 
implemented in the SWIFT package. We adopt values of Qp = 10^ and Q* = 10^° for our work. So tidal dissipation 
is dominated by that inside the planet. 



It is larger than the apsidal motion constant in binary studies 



(also written as k2) by a factor of 2. 



